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ABSTRACT 

Bulges in spiral galaxies have been supposed to be classified into two types: classi- 
cal bulges or pseudobulges. Classical bulges are thought to form by galactic merger 
with bursty star formation, whereas pseudobulges are suggested to form by secular 
evol ution due to spiral arm s and a barred structure funneling gas into the galactic cen- 
tre. iNoguchil (Il998l Il999l) suggested another bulge formation scenario, 'clump-origin 
bulge'. He demonstrated using a numerical simulation that a galactic disc suffers dy- 
namical instability to form clumpy structures in the early stage of disc formation since 
the premature disc is expected to be highly gas-rich, then the clumps are sucked into 
the galactic centre by dynamical friction and merge into a single bulge at the centre. 
This bulge formation scenario, which is expected to happen only at the high-redshift, 
is different from the galactic merger and the secular evolution. Therefore, clump- 
origin bulges may have their own unique properties. We perform a high-resolution 
iV-body/smoothed particle hydrodynamics (SPH) simulation for the formation of the 
clump-origin bulge in an isolated galaxy model and study dynamical and chemical 
properties of the clump-origin bulge. We find that the clump-origin bulge resembles 
pseudobulges in dynamical properties, a nearly exponential surface density profile, a 
barred boxy shape and a significant rotation. We also find that this bulge consists 
of old and metal-rich stars, displaying resemblance to classical bulges. These natures, 
old metal-rich population but pseudobulge-like structures, mean that the clump-origin 
bulge can not be simply classified into classical bulges nor pseudobulges. From these 
results, we discuss simila rities of the clump-orig in bulge to the Milky Way bulge. 
Combined with a result of lElmegreen et ahl (|2008l) . this pseudobulge-like clump-origin 
bulge could be inferred to form in clump clusters with a relatively low surface density. 

Key words: methods: numerical - galaxies: formation - galaxies: bulges - Galaxy: 
bulge - Galaxy: disc - Galaxy: formation. 



1 INTRODUCTION 



iKormendv fc Kennicutl (|2004l ') has suggested that bulges 
in spiral galaxies can be classified into classical bulges 
or pseudobulges. Classical bulges are thought to form 
through violent dynamical rel axation by galactic merger 
like elliptical galaxies ( e.g. iNakasato fc Nomotol 



Kormendv fc KennicuttI |2004|: ~ Naab fc Truiilld 



2002 



2ooe 



Rahimi et all I2OI0I : iHopkins et all bOld ). In other words 



both of classical bulges and ellipticals could be merger 
remnants of galaxies, therefore supposed to be struc- 
turally identical. Ind eed, their surface d ensity profiles 
follow the i?^/" law (|de Vaucouleurd Il948l '). their shapes 
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are rounder than pseudobulges and they have little net 
rotation. Since the galactic merger is expected to consume 
most of gas in bursty star formation, the systems evolve 
passively and tend to be a single population system of 
old stars. The rapid star formation leads the enhance- 
ment of Q-elements and m akes the bulge metal-rich 
(|Kormendv fc KennicuttI |2004| . and references therein) . It 
is well known that classical bulges and ellipticals compose 
an identical fundament al plane and luminosit y-mctallicity 
relation (e.g. iKormcndv fc lUingworth 19831: iKormendvl 
ll985l : |Jablonka et al.lll996l : lBurstein et al 
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Pseudobulges are discussed to form through secu- 
lar evolution caused by non-axisymmetric structures in 
a galactic disc, such as spira l arm s and a barred struc- 
ture (e.g. ICombes fc Sande"r3 Il98ll : IPfenniger fc NormanI 
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1990l: ICombes fc ElmegreenI Il993l: iBekki fc FreemanI [20o3 : 



Debattista et all |2004| : lAthanassoulal 120051 )7" These non- 
axisymmetric structures rearrange angular momentum of 
the gas in a disc, funnel the gas into the galactic cen- 
tre, which is presumed to arouse slow but long-lasting star 
formation resulting in a pseudobulge formation. Therefore, 
pseudob ulges display a wide variety in their stellar popula- 
tion. As lKormendv fc KennicuttI (|2004l ) has speculated that 
pseudobulges retain a memory of their discy origin, observa- 
tions indeed show disc-like natures of pseudobulges, such as 
an exponential surface density profile, a significant rotation, 
a flatter shape than those of classical bulges and frequent 
existence of nuclear structures (e.g. bar, spiral and ring). 

However, some bulges can not be simply classified 
into the two types above. The bulge of the Milky Way 
(MW) would be a typical example. Some observations 
have proposed a nearly exponential surface density pro- 
file, an oblate peanut shape (X-shape) and a significant 
rotation in the MW bulge, ind i cating the pseudobulge 
signa t ures (e.g. Kent et al.l 199ll: Kent 19921: Dwek et al] 



1 19951 : IWvse et all Il997l; [ Kormcndy fc Kennicut tI 12001^ 
and ref erences therein; lBabusia.ux ct al., .2010; Nataf et al.l 
2010l: [McWiUiam fc Zoccalil |201G|: 



Saito et al.l l201ll : 



Nataf fc Udalskil l201ll : iGonzalez etal] l2011bl )~The MW 



bulge also has characteristic properties of classical bulges. 
It has long been known that the MW bulge consists of old 
stars of which ages are suggested to be approxi mately 11 
- 13 G yr although uncert ai nty is still large (e.g. iTerndrud 
1 19881 : IWvse et all Il997l : iKormendy fc K ennicut tI |2004 
and reference s therein; ISahu et al.l 1200 6: Bcn sby et al.l 
l2010bl . l2011bl ). Additionally, rapid star formation episode 
were also suggested by metal-rich nature a nd an over- 
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ZoccaU et al.l 20061. 200a Rich et al.l 2007; Fulbright et al. 
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2010b 


. 2011bl; Johnson et al. 2011; Hill et al.l 2011^. Thus. 



these observations seem to indicate that the MW bulge 
can not be simply classified into classical nor pseudob- 
ulges. In addition, recent observations argued that some 
bulges, which were classified into pseudobulges, display 
exceptionally inactive star formation like classical bulges 
llDrorv fc Fished l2007l : iFisher et all |2009| ; iFisher fc Drorvl 
I2OI0I ). These inactive pseudobulges are also unclassifiable 
bulges. 

iNoguchil (|l998l . I1999I ) has proposed another bulge for- 
mation scenario; 'dump-origin bulge'. He demonstrated 
with numerical simulations that since galactic discs are 
expected to be highly gas-rich in early stage of the 
disc formation, clumpy structures form due to instabil- 
ity of the gas in the gaseous disc, which could also ex- 
plain some cl umpy gal axies observ ed in the high -redshift 



Universe (e.g. iNog uch? '1998", 'l999': 'imm oli et all l2004allbl ; 



Bournaud et al.|[2007uiAgcrtz et al...200a : .Aumer et al.ll2010l ; 



Ceverino et alll2010l . I2OI1I ). These galaxies are referred to 
as clump clusters (chain gal axies) i n the face-on (edge-on) 
view (e.g. Ivan den Bergh et al.,,1996 ; felmegreen et al.ll2004l . 
I2OO9I ; iGenzel et al.ll201ll ). iNoguchil ((l998l . Il999h suggested 
that these clumpy stellar structures fall into the galactic 
centre by dynamical friction and merge into a single bulge 
at the galactic centre, being clump-origin bulge. 



Clump-origin bulges form through 'mergers of the 
clumps' in a galactic disc, neither the galactic mergers nor 
the secular evolution. Therefore, properties of clump-origin 
bulges could be different fro m those o f the c onven tional 
ones. Since the simulations of INoguchil (| 19981 . I1999I ) have 
been performed with a low resolution, he could not discuss 
the detailed properties of the clump-origin bulge (Noguchi, 
private communication). A few foll owing studies discussed 
properties of the clump-origin bulge. lElmegreen et al.l (|2008l ) 
argued with numerical simulations that the clump-origin 
bulges should be classified into classical bulges, since the 
bulges in their simulations have a density profile near to 
the R^^* law, a thick shape, and a slow rotation in stellar 
component although the clumps are fast rotating before the 
bulge formation. They suggested that initial spin angular 
momentum of the clumps is lost to disc and halo when the 
clump-origin bulge forms. Recently ICeverino et al.l (|201ll ) 
studied dynamical state of gas in clumps migrating toward 
the galactic centre and demonstrated with both analytical 
approach and cosmological simulations that the gas com- 
ponent in most of clumps are highly rotation supported. 
Moreover, they discussed from their simulation results that 
more massive clumps are more highly rotation supported in 
the gas component and mentioned that clumps which have 
experienced a merger with other clumps tend to be highly 
rotationally supporte d ones (see their §4.2). These results of 
ICeverino et al.l (|201ll ) imply that the clump merger spins up 
the gas in the remnant clump. Since the clump-origin bulge 
is the final remnant of the clump mergers, we expect that 
the clump-origin bulges might be rotating fast like pseudob- 
ulges. However, they did not discuss the final structure of 
the bulge or dynamical state of the stellar component. 

We perform a similar (but more sophisti cated and much 
highe r resolution) numerical simulation to INoguchil (|l998l . 
Il999t ) using an isolated galaxy model and study the naive 
natures of clump-origin bulges in details. We reconsider the 
rotation of the clump-origin bulge mentioned above, finally 
discuss similarity of the clump-origin bulge to the MW 
bulge. We describe our simulation settings in §2 and the 
results of the bulge properties in §3, including dynamical 
state, stellar age, metallicity and star formation activity. We 
present discussion in §4 and summary of our results in §5 



2 SIMULATION 

The detailed description of our com puting schemes and in i- 
tial conditions has been presented in llnoue fc Saitoh! (|201lD . 
We describe it briefly here. 

We us e an A/'-body/SPH code, ASURA (jSaitoh et all 
|2008| . 120091'). Gas dv i iamic s is handled by the standard SPH 
scheme (|Monaghanl Il992l : ISpringell |201G| '). We here adopt 
an individual time-step scheme, the time-step limiter of 
ISaitoh fc Makind (|2009D to solve the evolution of the gas 
in the shocked regions correctly and the FAST method 
(|Saitoh fc Makind I2OI0I ) which accelerates the simulations 
of self-gravitating gas by integrating the gravitational and 
hydrodynamical parts with different time-steps. Gravity is 
solved b y a parallel t ree with GRAPE (GRAvity PipE) 
method (|Makind |2004| ) . We employ a barycentre approxi- 
mation and the tolerance parameter of 0.5. 

Cooling function of gas is assumed to be an optically 
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thin radiative cooling, depending on me tallicity and cover s 
a wide temperature range of 10 - 10* K (|Wada et al.ll2009l l. 
Feedback from a uniform far-ultraviolet radiation and type- 
II supernovae (SNe) are also taken into account. A gas par- 
ticle spawns a stellar particle of which mass is set to 1/3 
of the mass of the original gas particle under the criteria: 
1) pgas > 100 atm cm"^, 2) temperature T < 100 K, 3) 
V-'y< and 4) there is no SN around the particle. The lo- 
cal star formation rate is assumed to follow the Schmidt 
law (|Schmidtll 19591 ). The non-dimensional star formation ef- 
ficiency is set to C* = 0.033. Note that the global properties 
of star formation are fairly insensitive to the adopted value 
of Ct when one adopts a high mass resolution (~ IQ^Mq) 
and a high threshold density for star f ormation such as 
Pgas > 100 atm cm"^ JSaitoh et al.l 120081'). "W e assume the 
Salpeter initial mass function ( SalpeterTl955l ) ranging from 
0.1 Mq to 100 Mq on a stellar particle. We assume that 
stars heavier than 8 M© cause type-II SNe and thermally 
heat up ambient gas (see also §3.3). 

Our initial c ondit i on fo llows the spherical model of 
iKaufmann et al" I (|2006l . 120071 ) that was used to study the 
formation of disc galaxies in an isolated environment. We as- 
sume a n equilibrium system with the Navarro-Frenk- White 
profile (jNavarro et al. 1 19971 ) with a virial mass Mvir = 
5.0 X 10^^ Mq, a virial radius rvir = 1-67 x 10^ kpc and 
a concentration parameter, c = 6.0. Baryon mass fraction 
of the system is set to 0.06. Initially all of baryon is pri- 
mordial (zero-metal) gas having a virial temperature. The 
density profiles of dark matter (DM) and gas components 
are initially the same. The gas sphere has a rotation fol- 
lowing a specific ang ular momentum distribution of j oc r 
iBuUock et al.l l200ll ) normalized by a spin parameter of 
A = 0.1. The DM halo is represented by 10^ particles with 
a softening length of edm = 8 pc, whereas the gas in the 
halo is expressed by 5.0 x 10^ particles with tgas ~ 2 pc. 
The masses of DM and gas particles are 4.7 X lO'^M© and 
6.0 X 10^ M0, respectively. 

Although we solve the wide temperature range of the 
inter-stellar matter (ISM), T = 10 - 10* K, we adopt nei- 
ther temperature nor pressure floor. This is because we think 
these artificial suppression techniques of the low Jeans mass 
objects are unnecessary. It is widely believed that insufficient 
mass resolution may lead to unphysical fragmentation in 
the ISM. To avoid this 'unphysical fragmentation', one may 
adopt temperature/pressure floors which can maintain the 
condition of Mjeans > Ai'nb X mgas, where Afjeans is the Jeans 
mass, Nnb is the number of neighbouring SPH particles and 
rrigas is the mass of a gas particle. Since the original claim of 
the unphysical fragmentation was do ne based on their simu - 
lations during non-linear evolutions (|Bate fc Burkert|[l997l ). 
it is hard to clarify what was the true origin of the unphys- 
ical fragments. This p roblem has been first pointed out by 
iBate fc Burkerdll 19971 ). Since their evaluation was based on 
their simulation during non-linear evolutions, it is difficult 
to clarify what was the actual origin of the unphysical frag- 
ments. 

iHubber et al] (|2006l ) showed with well-conceived nu- 
merical tests for this problem that the growth of a linear- 
regime perturbation in a self-gravitating fluid is just slower 
than that of the analytical solution when the mass resolu- 
tion is insufficient to express the local Jeans mass (Mjcans > 
A^/rm7i6 XTTigas). They did not find any unphysical fragmenta- 



tion in simulations with insufficient mass resolution. There- 
fore, there is no need to worry about the unphysical frag- 
mentation in simulations with insufficient mass resolution 
without temperature/pressure floor. Hence, the result that 
the disk holds a number of small clumps especially in early 
phase in our simulation is not a numerical artifact (see, §3). 

The effect of the insufficient mass resolution in our sim- 
ulation would appear as the suppression of growth rates of 
clump seeds. However, this is not crucial for our results be- 
cause (1) the dynamical time of seeds we may drop in our 
simulation is short because they are compact (2) the final 
mass of the clumps are much larger than the resolution limit. 



3 RESULTS 

We run the simulation for 6 Gyr and display snapshots of 
stellar and gas distributions and star forming region in Fig. 
[1] The simulated galaxy settles into a stable state after 5 
Gyr. 

At t < 2 Gyr, many clumps form in the disc onto which 
the surrounding gas is falling. Active star formation can 
be seen only in the clumps, whereas little star formation 
is found in the inter-clump (disc) region (the third row). 
This clumpy phase is seen in previous numerical studies and 
suggested to explain the clump clusters and cha in galax- 



ies o b serve d in the high-redshif t Universe (e.g. Noguchil 
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In our 


simulation model, since the system is lacking 



continuous gas supply supposed in the cosmological context, 
the gas accretion ceases at f ~ 2 Gyr. After the cessation 
of gas supply, only a small number of tiny clumps form in 
the disc, which are fragile to tidal force and these disrupted 
clumps build up a stellar disc. The other massive clumps can 
survive, merge with one another and grow to more massive 
ones. In this phase, a central density slope of the DM halo 
is turned co red by the reaction of the dynamical friction on 
the clumps (|lnoue fc Saitohll201ll ). The clumps are sucked 
into the galactic centre by dynamical friction and finally 
m erge into a sing l e 'clu mp- origin' bulge as first suggested 
bv lNoguchil (|l998l . ll999l ). At the final state (t = Q Gyr), 80 
per cent of the initial gas had converted into stars; the total 
stellar and gas masses are 2.5 x 10^° M© and 5.0 x 10® Mq, 
respectively. There is little gas in the halo region. The central 
density cusp of the DM halo is revi ved due to the gravita - 
tion al potent i al of t he massive bulge (|lnoue fc Saitohil201ll ). 

iNoguchil 1 19961 ) has first suggested a thick disc forma- 
tion scenario by dynamical h e ating on a disc by clumps, 
and lElmegreen fc ElmegreenI (|2005h observationally dis- 
cussed the exis t ence of a thick disc in chain galaxies. 
iBournaud et all (|2009l ) demonstrated with numerical sim- 
ulations that a disc forming in clump clusters is a thick 
disc, whereas a thin disc forms via gas accretion at the 
latter epoch than that of the thick disc formation. Since 
our simulation model does not take such a late stage ac- 
cretion into account, we indeed find in our simulation that 
the disc is dominated by the thick disc component with a 
scale height of 1.41 kpc and a scale radius of 9.85 kpc, while 
a thin disc hardly forms. Therefore, in this paper, we con- 
sider the disc in our simulation to be a single-component 
thick disc which forms via the clumps. However, it does 
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Figure 1. Stellar and gas density maps in tiie central 60 X 60 kpc region. Tlic first and second rows indicate the surface density 
of all stars from the edge-on and face-on views, respectively. The third row plots the density of stars younger than 40 Myr as star 
forming region. The fourth row shows the volume density of gas in a slice of the disc plane. Time in unit of Gyr is indicated on 
right bottom in each panel. 



not mean to reject other possible scenarios of thick disc 
form ation, such as dissipativ e coUapse (e.g. [ Burkert et ahl 
Il992l ) , thin disc heati ng (e.g. Quinn et al.l 19931 ). accretion 
of dwarf galax ies (e.g. Abadi et al.l 1 20031 ') . muhiple gas-rich 
min or merger |Brook et al]l2 004. 200^) and radial migration 
(e.g. lSellwood fc Binnevll200a '). We will report aU the details 
of the thick disc formation and kinematics in our simulation 
in a forthcoming paper (Inoue & Saitoh, in prep). 

In Fig. [2l we show the star formation history of the 
entire galaxy. Since the gas is continuously falling onto the 
galaxy {t <2 Gyr), the system keeps star formation rate 
of ~ 10 Mq yr~^. After the gas infall phase, the galaxy 
asymptotically settles to the low star formation state with ~ 
1-0.1 Mq yr"'^. When the stellar clumps coalesce with one 
another, the star formation activity is temporarily activated 



in the clusters, resulting in appearance of spiky features in 
the late history of star formation in Fig.H 



3.1 Dynamical state 

3.1.1 Density profile 

In order to identify the bulge type, i.e. classical or pseudob- 
ulges, the Sersic profile fitting of a surface density profile of 
bulge stars is frequently used: 



E(i?) = Eoexp 



(i?o) 



(1) 



where Eq and 7?o are the central surface density and the 
scale radius of the bulge, respectively, and n is the Sersic 
index. If n = 1, the profile is consistent with an exponen- 
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Figure 2. Star formation rate (SFR) as a function of time. The 
entire region of the galaxy is considered. 
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Figure 3. Stellar surface density profile from the face-on view at 
the final state of the simulation (t = 6 Gyr). Minimum fitting 
by the Sersic profile with the index, n = 1.18, is also plotted. 
The centre of the system is defined to be the densest point in 
the stellar density distribution. The insert diagram illustrates the 
same density profiles with the linear abscissa. 



tial profile implying a pseudobulge. If w = 4, the profile 
follows the de Vaucouleurs' R}^'^ law (|de Vaucouleurslll948l ) 
implying a classical bulge. The bo undary is empirically set 
to n ~ 2 with little overlap (e.g. Kormendv fc KennicuttI 
l2004l : iDrorv fc Fisheill2007l : iFisher fc Drorvll2008l ). 

In Fig. [31 we plot the azimathally averaged surface 
density of stars at the final state and the result of the 
Sersic profile fitting. The fitting by minimum is given 
by Eo = 1.37 x 10* M0/pc^ Ro = 0.21 kpc, and the Sersic 
index, n = 1.18. This value of the index, n, means the bulge 
in our simulation has a nearly exponential surface density 
profile and indic ates a pseudobulge-like density structure. 
Contrary to this. lElmegreen et al.l (|2008l ) has indicated with 
numerical simulations that their clump-origin bulges have a 
Sersic index larger than 2 and the morphology of classical 
bulges. We will discuss this contradiction in §4.4. 

The Sersic profile can not fit the simulated profile out- 
side i? ^ 1 kpc. We find the transition radius between the 
disc and the bulge components at i? 1 kpc on the sur- 
face density. Thus, we define the bulge region to be inside 



R — 1.0 kpc. The bulge has a mass of 4.5 x 10 M© and a 
half-mass radius of 450 pc. Ratio of the bulge mass to the 
total stellar mass, B/T, is 0.1^. However, it must be kept in 
mind that the galaxy in our simulation does not have a thin 
disc. Since a thin disc accounts for not a small fraction of 
stellar mass of a re al galaxy, the B/T mus t be overestimated 
in our simulation. iFisher fc Drorvl (|2008l ) has observed the 
mean value of B/T to be (B/T) = 0.16 ± 0.05 among pseu- 
dobulges and {B/T) — 0.41 ±0.11 among classical bulges. 
Therefore, the clump-origin bulge in our simulation is near 
to pseudobulges on the value of B/T. 

It should be noted that the B/T could be sensi- 
tive to simulation setup. We speculate that the B/T 
would become smaller if the gas accretion onto the disc 
is slower, because wavelength at which insta bility first ap- 
pears becomes smaller in a gas-poorer disc (|Toomre|[l963 : 
iBinnv fc Tremain^ l2008l l. which allows to form smaller 
clumps preferentially. They tend to be disrupted by the 
galactic tide and/or lead to result in a smaller clump-origin 
bulge. On the other hand, if we assume a smaller spin pa- 
rameter in the initial condition, the gas forms a smaller disc 
and the gas is mor e condensed in the disc than that in the 
current simulation (|Fall fc Efstathioulll98Gl : lBaron fc White! 

In such cases, the wavelength of the instability be- 
comes larger, the larger clumps would form and result in 
the formation of a larger bulge. In the dense and small disc, 
since the timescale of the clump migration becomes much 
shorter, the bulge formation would be finished earlier than 
the current simulation. 



3.1.2 Shape 

In Fig. 21 we show stellar surface density maps of the bulge at 
the final state in our simulation. As seen in the image from 
the face-on view (the top panel), the bulge has a bar- like 
elongated feature. This 'nuclear ba r' in a bulge is a common 
featu re among pseudobulges (e.g. iKormendv fc KennicuttI 
l2004h . EUipticitjE] of the bulge measures 0.55 from the face- 
on view. 

As seen from the edge-on view (the bottom panel), it 
clearly appears that this bul ge is a boxy bulge, indicating 
a pseu dobulge signature (e.g. iKormendv fc Kennicuttll200"^ 
lAthana ssoula 2005,). EUipticity of the bulge measures 0.67 
from the edge-on view. 



3.1.3 Rotation 

Pseudobulges are generally observed to ha ve a significant 
rotation (e.g. Kormendv fc KennicuttI [2004I . and references 
therein; I Williams et al.ll201lh . On the other hand, classi- 
cal bulges are generally presumed not to rotate signifi- 
cantly although some classical bulges show a net rot ation 
(|Kormendv fc Illingwor"thlll982l : ICappellari et al.ll2007h . 

In Fig. [5l we illustrate a map of the mean line-of-sight 
(LoS) velocity of stars from the edge-on view (perpendicular 



^ For example, the MW has a v alue of B/T = 0.19 ± 0.02 
llKent et al.lll99ll : iDwek et al]|l995^ . 

^ In this paper, the ellipticity is defined to be 1 — 6/a of the 
isodensity contour at which the surface density is equal to So/2, 
where a and b arc lengths of major- and minor-axes. 
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Figure 4. Stellar surface density maps in central 3x3 kpc region 
at the end-state of the simulation (t = 6 Gyr). The top panel is 
from the face-on view and the bottom one is from the edge-on 
view. The horizontal axes (rr-axes) coincide with the major-axis 
of the bulge in both panels. The contour levels in the images arc 
chosen to highlight the bulge shape. 



to the major-axis) at the final state. The value is normal- 
ized by the LoS velocity dispersion interior to the half mass 
radius of the bulge (450 pc), ao = 81.0 km s~^. This bulge 
shows a significant rotation. The value of Vniax/o"o — 0.9 
of this bulge means that the rotation (spin) is not negligi- 
ble but comparable to the velocity dispersion, where Knax 
is the maximum mean LoS veloci ty in the bulge. This is 
also a pseudobulge-like signature (|Kormendv fc KennicuttI 
I2OO4I . and references therein). The direction of the rotation 
coincides with the disc rotation. 




<V, 



Figure 5. Moan LoS velocity map of the stellar component in the 
central 3x3 kpc region at the final state of the simulation (t = 
6 Gyr) from the angle edge-on and perpendicular to the major- 
axis of the bulge. The value is normalized by the LoS velocity 
dispersion interior to the half mass radius of the bulge, ag = 
81.0 km s^i. 
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Figure 6. Rotation velocity curves at different heights above the 
disc plane, z. The abscissa coincides with the major axis of the 
bulge. 



Recently, observations of some pseudobulges have 
shown that rotation al velocities do not dep end on a height 
within the bulges (jWilliams et al.l 



iepi 



20111 . and references 



Howard et alll2009l ). 



therein), including the MW bulge 
'cylindrical rotation'. However, regardless of our results that 
the clump-origin bulge follows the pseudobulge-like natures 
and the boxy morphology, we do not find a cylindrical rota- 
tion in the clump-origin bulge simulated in this work. Fig. [6] 
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Figure 7. Averaged formation time of star particles as a function 
of tlie galactocentric distance. The values are azimathally aver- 
aged and a dispersion, a, is calculated in each radial bin. This 
is a formation time of stars, not age. We define the start of the 
simulation to be the zero-point of the formation time. 



shows that the rotation velocity curve s tend to be flatter at a 
higher latitude of the bulge. However . I Williains" et al.l(|201ll ) 
has observationally shown that boxy/peanut-shaped bulges 
do not necessarily involve a cylindrical rotation. Therefore, 
their result implies that the cylindrical rotation may not be 
a congenital nature of the boxy /peanut-shaped bulges. 

Recently, ISaha et al.l (|201ll ) indicated with a numerical 
simulation that interaction with a barred structure in a disc 
can spin up a non-rotating bulge and the bulge acquires a 
boxy/peanut shape and a cylindrical rotation. If this mecha- 
nism is efficient, the cylindrical rotation can be built up after 
the bulge established and may not be a strong indicator for 
the bulge formation mechanism. 



3.2 Age 

As we noted in §1, classical bulges are generally an old struc- 
ture because the classical bulges form by past galaxy merger 
events. Therefore, their stellar evolution tends to be passive 
after the last merger. As for pseudobulges, their star for- 
mation is generally still going on (although see also §3.4) 
because non-axisymmetric structures continue to funnel gas 
into the central bulge. Hence, stellar age in a bulge is also 
another important indicator to identify the bulge type. 

In Fig. [T] we plot the mean and dispersion of the for- 
mation time of stars in our simulation. As seen in the figure, 
the bulge [R < 1 kpc) is younger than the thick disc only 
by ~ 1 Gyr. It is because the clumps in the clump cluster 
retain star formation activity inside them (the third row of 
Fig. m and migrate to the galactic centre while forming new 
stars. The formation time of disc stars are almost indepen- 
dent from the distance from the galactic centre. 

However, the difference in the formation time between 
the bulge and the disc, ~ 1 Gyr, is so small in comparison 
with the age of the Universe that we could presume the 
bulge to form simultaneously with the disc. As mentioned 
above, the disc in our simulation is a thick disc, does not 
have a thin disc. Moreover, this bulge formation scenario 




Figure 8. Radial mctallicity distribution at t = 6 Gyr. The 
dashed lines indicate the range of itlo" dispersion. The value of 
solar mctallicity, Zq, is set to 0.02. 



is expected to happen only in the high-redshift Universe. 
Therefore, we could consider that this clump-origin bulge is 
an old structure, as old as the thick disc. 

To be an old structure does conflict with the morpholog- 
ical and dynamical similarities of this clump-origin bulge to 
pseudobulges mentioned in §3.1. Therefore, we suggest that 
the clump-origin bulge dynamically resembles pseudobulges 
but consists of the old stars like classical bulges. Hence, the 
clump-origin bulge can not be simply classified into classi- 
cal bulges nor pseudobulges, but shows peculiar properties. 
However, we have to note that the bulge may be rejuvenated 
by the late accretion to form a thin disc. 



3.3 Metallicity 

Before discussing metallicity, we should describe here how 
we treat chemical evolution in our simulation. Our simula- 
tion code tracks only the total metal abundance, Z, and does 
not follow the evolution of each element. The initially pri- 
mordial gas is contaminated with heavy elements by type-II 
SNe. We did not take type-la (or other) SNe or stellar wind 
from intermediate mass stars into account in our simulation. 

We show a radial metallicity distribution in our simu- 
lation in Fig. [S] Stars of the clump-origin bulge are more 
metal-rich by 0.2 dex than stars of the inner region of the 
thick disc. Such a high metal abundance in the bulge is a 
natural outcome of the formation of clumps which form from 
a collapsing gas cloud in the highly gas-rich disc and the vi- 
olent star formation in the clump mergers. 

Recent observations have shown that there is no or lit- 
tle chemical distinction between the MW bulge and the 
thick disc stars in [a/Fe] vs. [Fe/H] distributions despite 
exhibiting clear differences in their metallicity distribu- 
tion f unc tions jMelendez et all l2008l: iBensbv et all |2009| 
2010al |h l. l2011bl: lAlves-Brito et all l20ld : iGonzalez et al] 
2011al: I Hill et al] l201lh . However, observation results of 
Lecureur et all (|2007l ') and iFulbright et all (|2007l 'l indicated 
clear separations of [a/Fe] between the thick disc and bulge 
stars in the MW. Observations of the chemical distinction 
between the MW bulge and the thick disc stars are still 
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Figure 9. Vertical mctallicity distribution at tlie bulge regions 
(la;|=0.0 kpc and 0.5 kpc) and disc regions (|x[=3.0 kpc and 
10.0 kpc). Values are averaged along LoS. Thin lines (red solid 
and green long-dashed) indicate the metallicity distribution of the 
bulge stars only {R < 1 kpc). 



controversial because of only a handful number of reliable 
sample stars and lack of accuracy of distances and proper 
motions to distinguish the bulge and disc stari|f|. Our simu- 
lation predicts the higher metallicity in the bulge than the 
thick disc. More convincing observations are highly antici- 
pated to test such a scenario. Note, however, that there are 
ma ny ways to bu il d up t hick discs, as we noted in §2. 

iRuchti et al] (|201lf ) recently observed that the MW 
thick disc shows a small radial gradie nts in [Fe / H] ra - 
tio, -1-0.01 ±0.04 dex kpc~\ [All cndc P rieto et all (|2006l ). 
Ilvezic et all (|2008l ) and iBensbv et al.l (|2oTlal ) also observa- 
tionally confirmed no or little radial [Fe/H] gradient in the 
thick disc stars. Our result in Fig. [8] shows a small neg- 
ative gradient, —0.014 dex kpc^^, in metallicity, Z, from 
R = 2.0 kpc to 15.0 kpc, which is not inconsistent with 
these observations in the MW within the error. 

It has been indicated that there is a clear vari- 
ation i n metallicity along the minor axis of th e MW 



bulge (Ivan d en Bcrgh & Herbst Il974l : iTcrndrud 
Frogel et al1]l990: Tcrndrup et al.lll990l: [Zoccali ct all 



1988 



2008 



Johnson et al.l 201ll : iGonzalez et al.l l2011bl ). IZoccali et al.l 



(|2008l ) observed the vertical gradient of the mean [Fe/H] 
ratio to be —0.6 dex kpc^^. We here show metallicity distri- 
butions along the vertical direction in our simulation in Fig. 
[2] This figure clearly indicates a vertical metallicity gradient 
along the minor axis in the clump-origin bulge, z < 1.0 kpc 
of the solid red line. Contamination by the disc stars is neg- 
ligible in the bulge region (the thin lines). In our result, 
the metallicity keeps increasing eve n in the central r egion 
of the clump-origin bulge although I Rich et al] (|2007l ) has 
observed no evidence for the vertical gradients of [Fe/H] 
and [of/Fe] in the central region of the MW bulge, between 
{l,b) = (l°,-4°) and {l,b) = (0°,-l°) in the Galactic coor- 
dinate. 

In the MW thick disc, iRuchti et al] (|201ll ) observed 



Figure 10. Specific star formation rate at the final state of the 
simulation (t = 6 Gyr). Filled and open symbols indicate values 
averaged in the last 0.5 Gyr and 40 Myr, respectively. Circles 
and diamonds indicate the values of the bulge and disc region, 
respectively. The open symbols are missing data points at inner 
disc region, 1 kpc < R, < 6 kpc, since no star particle was born 
in this radial range in the last 40 Myr. 



that the vertical gradient of [Fe/H] ratio is very small, 
— 0.09±0.05 dex kpc^^. In our simulation, the vertical gradi- 
ent of Z is small, ~ —0.05 dex kpc~^ in the thick disc region. 
Additionally, we find that the vertical metallicity gradients 
are almost independent of the galactocentric distance. This 
is represented in Fig. |§1 The blue dotted line (inner disc re- 
gion) and the purple dot-dashed line (outer disc region) are 
almost parallel. 



3.4 Star formation activity 

Classical bulges are supposed to have used up their gas 
and/or have a central massi ve black hole to hea t up and 
prevent gas from accreting (|Okamoto et al.l |2008| ) . There- 
fore, they are generally inactive in star formation. Pseudob- 
ulges are still forming stars actively, which is observation- 
ally verified by star forming nuclear rings nuclear spirals, 
bars and ovals, implying sim ilarities to disc galaxi es (e.g. 



Kormendv fc Kennicutt 20041 . and references therein: iFished 
20061 : iFisher fc Drorvll2010l ). 

In Fig. 1101 we show the specific star formation rat^ in 
the bulge and the thick disc of our simulation. The clump- 
origin bulge is little forming stars at the final state. This 
bulge has very quiet and weak star formation activity. In the 
thick disc, the specific star formation rate becomes larger as 
R increases. The innermost region of the disc is less active 
than the bulge region. There is little difference between the 
values averaged in the last 40 Myr (the open symbols) and 
0.5 Gyr (filled symbols). Thus, the inactive star formation 
is not transient. This halted star formation in the bulge is a 
typical feature of classical bulges. 

Some bulges classified into pseudobulges are discov- 
ered to show inactive star formation like classical bulges 



3 iHill et al.l feOllh argued that using different reference stars in 
calibration was the cause of the confiict. 



^ The specific star formation rate is defined in this paper as a 
star formation rate in a radial bin divided by the mass of the bin. 
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JProrv fc Fished l2007l : iFisher et all |2009| : iFisher fc Drorvl 
These bulges have a smaU value of the Sersic in- 
dex, n < 2, despite the fact that their star formation activ- 
ity is as weak as classical bulges . These are called ina ctive 
pseudobulges' l|Fisher et al.ll2009l : [Fisher fc Drorvll2010l '). Al- 
though the bulge in our simulation is less active than the 
observed inactive pseudobulges by ~ 1 dex in the specific 
star formation rate (~ 10~^^ yr~^ in the observations), the 
clump-origin bulge in our simulation resembles to the in- 
active pseudobulges: the nearly exponential surface density 
and the inactive star formation. 

We have to note that it could not be fair to directly 
compare the bulge in our simulation with the inactive pseu- 
dobulges in observations. The inactivity of these observed 
pseudobulges may be due to feedback from an active galactic 
nucleus, while our simulation does not take it into account. 
Moreover, the simulated galaxy is lacking the later gas ac- 
cretion which would play an important role for the thin disc 
formation. The star formation activity in the clump-origin 
bulge would be affected by the gas accretion: a dynamically 
cold thin disc forms non-axisymmetric structures such as 
spiral arms and a galactic bar, which can supply the bulge 
with gas more or less and would activate the star formation 
in the bulge. 



4 DISCUSSION 

4.1 Validity of the model settings and comparison 
with observations 

We set the larger spin parameter than t he mean value of cos- 
mological simulati ons, A ~ 0.035, (e.g. iBuUock et al]|200ll : 
iKimm et al1l201ll ) although observed values among high- 
redshift galaxies are stil l controversial (e.g. iBouche et al.l 
I2OO7I : iDutton et alllioTl] ). As we mentioned in §3.1.1, the 
large spin parameter would lead the formation of a large disc, 
small clumps and a small bulge. lElmeereen fc ElmegreerJ 
(|2005l ) observed that diameters of galactic discs range 13 
- 29 kpc among 10 clump cluster samples of which masses 
are 0.5 - 25 xlO^^M©. The numbers of clumps in the sam- 
ple are 5-14 per galaxy and the averaged clump masses 
in the clump cluster range 0.21 - 1.36 xIO'^Mq. The disc 
in our simulation is larger than their observational result 
by a factor of 2 or 3 and the surface density of inter- 
clump region is lower than the observed values (see also 
lElmegreen et al.ll2009l ). Therefore, our simulation may not 
be a well-reproduced model of the high-redshift galaxies in 
the current observations on the points of the disc size and the 
surface density, although our simu lation does not seem in- 
consis tent with the observation of lElmegreen fc Elmegre"enl 
(|2005l) on the points of the number and mass of the clumps 
(compare with 2 — 3 Gyr in Fig[T]). In our simulation, the 
mass resolution of gas, which corresponds to the minimum 
resolved Jeans mass, is approximately N^b x JTigas — 2 x 
10^ M0. T his is much smaller than the me an clump masses 
observed in lElmegreen fc Elmcgrccn' ('2005') and comparable 
to the smallest clumps in .Elmcgrccn ct al. (2009). Therefore, 
the mass and number of clumps we obtained in our simula- 
tion are do not suffer from the resolution effects. 

However, current high-redshift observations still have 
difficulties of limited angular resolution and detectability. 



which would have caused underestimates of the number 
of clumps and the disc diameter. Moreover, low luminos- 
ity (dum£_clu^tera might have been missed. As discussed 
in lElmegreen et al.l (|2009l ). the current high-redshift obser- 
vations would be biased to bri ght parts in bright galaxies. 
lElmegreen fc ElmegreerJ (j2005f l discussed that their clump 
cluster samples have a too high inter-clump surface density 
to be a typical nearby spiral galaxy like the MW. Although 
clump clusters with a low surface brightness would not be 
the most representative cases for the high-redshift clump 
clusters, it may be possible that the low surface brightness 
clump clusters had been present but missed in the observa- 
tions, and may be a progenitor of the nearby spiral galaxies. 
Our simulation result may corresponds to such a low lumi- 
nosity clump cluster, which suggests the old and metal-rich 
pseudobulge formation from the clumps. 

The longevity of the clumps in our simulation should be 
also examined carefully. In our simulation, the clumps take 
~ 5 Gyr to merge into the bulge. The clump migrations 
in the previous simul ations have accompli s hed in a shorter 
time, ~ 1 G yr (e.g. Ilmmeli et al.ll2004al: iBournaud et al.l 



2OOI I2OO9I: Elmegreen et al.l 120081 : lAumer et all bOloT 



Ceverino et al.l2011 ). The main causes of this difference are 



the disc size and the clamp mass. The clumps in our simu- 
lation have to migrate a long distance to the galactic centre 
because of the larger disc scale length than that in the pre- 
vious works. In addition, the clump masses seem somewhat 
smaller than these in the previous works. Therefore, the dy- 
namical friction on the clumps is weak in our simulation. 

In such a case of the slow evolution with the long-lived 
clumps, the clump clusters can not finish the formation of 
disc and bulge in the early Universe. If the evolution time 
scale is ~ 5 Gyr like our simulation, the clumps remain 
in galaxies until after the redshift z ~ 1. Observationally, 
there exists clum p clusters even after the redshift z ~ 1, and 
lElmegreen et al. (2 009 ') observed that the clump clusters in 
the lower-redshift Universe are less massive than these in 
the higher-redshift Universe although they discussed that 
this 'downsizing' is due to a detection limit in the obser- 
vation and a smaller sampling volume in the lower-redshift 
Universe. If this downsizing is real, our simulation may cor- 
respond to such less massive galaxies evolving slowly until 
the low-redshift. 

Di sruption of th e clum ps should be also discussed. Re- 
cently, iGenzel et al.l ()201ll ) has observed strong outflows 
from some clu mps, which a r e fast er than the escape velocity 
of the clump. iGenel et al] (j2010l ) demonstrated the clump 
disruption due to the superwinds implemented in their nu- 
merical simulations by a p arametrised wind model, although 
iKrumholz fc Dekej lj2010l ) had analytically discussed that 
such a strong feedback would not be realistic. The SN feed- 
back model adopted in our simulation can not repro duce 
such destructive superwinds. IKrumholz fc Dekell (|2010h sug- 
gested that the most probable engine of the superwinds is ra- 
diation pressure from newly formed massive stars. However, 
most of simulations for galaxy formation have not been im- 
plemented with the radiation pressure from the massive stars 
and o ur simulation did not either. The result of lGenel et al.l 
(|2010h implies that taking the radiative feedback into ac- 
count may drastically change the picture of the simulations 
for clump clusters. 
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4.2 Similarity to the MW bulge 

The MW bulge does not follow the properties of classical 
bulges nor pseudobulges. The MW bulge shows a nearly ex- 
ponential profile, an oblate peanut shape (X-shape) and a 
significant rotation, which are similar to pseudobulges. At 
the same time, it has old stars, a high metallicity and weak 
star formation, which are similar to classical b ulges (e.g. 
IWvse et all 1 19971 : iKormendv fc Kennicuttll2004 and refer- 
ences therein). These properties of the MW bulge are consis- 
tent with the clump-origin bulge obtained by our simulation 
although the absence of cylindrical rotation is inconsistent 
with the MW bulge. 



be 



Bulge stars of the 
a-element enhanced 



MW are observed to 



19941: IKormendv & Kennicutt 


I2OOJ. and references 


therein; Zoccali et al 
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2007: Melendez et all 20081: 
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: ImU et al.l 2011). Interestindv. 


Zoccali et al.l (|2006|) observationallv discussed that an 



[0/Fe] vs. [Fe/H] distribution of the MW bulge stars agrees 
very well with the result of the chem o-dynamical simulatio n 
for clump-origin bulge formation in llmmeli et ahl (l2004al ). 
On the other hand, contrary to IZoccali et al.l ), the 

recent observation of Baade's window by Hill et al. ( 201 ll ) 
showed considerable disagreement in a met allicity distri- 
bution function with the simulation result of llmmeli et al.l 
(|2004al ). These chemical abundances will provide further 
constraints on the bulge formation scenario. 

We assumed an idealized simulation model, a collaps- 
ing gas sphere in a Navarro-Frenk- White model halo, in or- 
der to study naive natures of a clump-origin bulge. As we 
noted repeatedly in previous sections, our model would need 
the l ater gas accretion to form a thin disc (|Bournaud et al.l 
I2OO9I ). Moreover, in the cosmological context, real galaxies 
suffer galaxy merger more or less. These physical processes 
following the clump-origin bulge formation could alter the 
naive natures of the clump-origin bulge. For example, as we 
noted in §3.1.3, a bar red structure may induce a cylindrical 
rotation of the bulge jSaha et al.ll20ll1 ). 

Recently, some observations showed bimodal metal- 
licity distributions and kinematics among MW bulge 
stars, implying a multiple formation process of the MW 



stars, implying a multiple tormation process 01 tne ivivv 
bulge (iBabusiaux et aL|]2010l: iBensby et all r2010bl. l2011bl: 



iHill et all I2OIII ). iBabusiaux et al.1 (|2010D and iHiU et all 
I 2OIII) suggested that the MW bulge is a complex bulge of 
a classical bulge and a pseudobulge. Thus, bulge formation 
in real galaxies (including the MW) might be complex and 
proceed through multiple scenarios. Some bulges would be 
a composite struct ure as proposed by some observations of 
extra-galaxies (e.g.lErwin et al.l | 2003l: Kormendv et al.ll2010l : 



IKormendv fc Barentinell2010l : iFiTher fc Drorvll201ll ) 



Such unclassifiable bulges like the MW bulge (old pseu- 
dobulge) are a lso obs erved in some oth er disc galaxies (e.g. 
iBica fc Alioinlil987, : ,Peletier et al.lll999l ). Our simulation re- 
sults suggest that clump-origin bulges could be a possible 
origin of these old pseudobulges. In addition, as we discussed 
in §3.4, the inactive pseudobulges may be also clump-origin 
bulges. 



4.3 Exponential bulges in the high-redshift 
Universe 

It is very difficult to determine whether a bulge observed 
in a distant disc galaxy is clu mp-origin or not. Although 
IKormendv fc KennicuttI (|2004l ) has argued that the secular 
evolution could build some pseudobulges quickly > 5 Gyr 
ago, we expect that the secular evolution would not be able 
to build up a pseudobulge in a premature disc in early phase 
of the galaxy formation. It is because the ^ secular^ evolution 
has to postdate accomplishment of their disc formation and 
takes a long time to form a pseudobulge. Therefore, we can 
expect that a galactic bulge observed in the high-redshift 
Universe would be either a classical bulge or a clump-origin 
bulge. As we showed in §3.1.1, the clump-origin bulge has 
a nearly exponential surface density profile, which is distin- 
guishable from the R^^* law profile of the classical bulges. 
Thus, if the density profile of the bulges in clump clusters 
could be observationally resolved, we would be able to iden- 
tify the clump-origin bulges by identifying their exponential 

den sity profile. 

lElmegreen et al.l (|2007l ) has observationally indicated 
that comoving spatial number density of clumpy galaxies 
(clump clusters and chain galaxies) in the high-redshift Uni- 
verse (1 < z < 4) are consistent with the density of spiral 
galaxies in the nearby Universe and suggested that quite a 
few galaxies of current spirals used to be a clump cluster or a 
chain galaxy. This implies that disc galaxies hosting a clump- 
origin bulge may not be rare. However, the idea that all disc 
galaxies must have formed a clump-origin bulge contradicts 
the observational fa ct that the number of bulge-less spiral 
galax ies is not small (IKormendv et al" ]|20ld:lFisher fc Droryl 
I2OIII ). As we noted in H.l. lElmegreen fc Elmegreen! (|2005l ) 
had discussed that the observed clump clusters could not be 
the nearby spiral galaxies because of the too high surface 
density. 



4.4 Comparison with previous studies 

So far, there have been only a few studies discussing 
the natures of c l ump-o rigin bulges. As we noted in §3.1, 
lElmegreen et al] (|2008l ) has argued with numerical simu- 
lations that clump-origin bulges have a similar dynamical 
state to classical bulges, such as a large Sersic index, n ^ 2 
and a slow rotation, V/ao ~ 0.4 — 0.5. Their results indicat- 
ing classical bulge-like signatures poses a discrepancy with 
our results. This disagreement may imply that the properties 
of clump-origin bulges depend on initial conditions of sim- 
ulations. Our simul ation model adopted th e collapsing gas 
sphere model, while lElmegreen et alJ (|2008l ) used a uniform 
density disc model as their initial condition. As we discussed, 
the size and number of the clumps would depend on the ini- 
tial spin parameter. The size of the clumps may aff'ect the 
re laxation process in the clump mergers. In the simulations 
of lElmegreen et al] (|2008h , fewer and larger clumps form in 
the clump cluster and the resultant B/T ratio is somewhat 
larger than in our simulation. 

Simulation method may also b e a possible cause of the 
disagreement. The simulations of lElmegreen et al] (|2008l ) 
were performed by a particle-mesh sticky-particle code 
which mimics hydrodynamics by inelastic collisions be- 
tween particles designated as a gas particle. On the other 
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hand, we describe the multi-phase interstellar medium 
with our SPH code which treats radiative cooling, far- 
ultraviolet heating and SN feedback with hydrodynami- 
cal equations of motion. Since the sticky particle method 
is only an imitation of hydrodynamics, the evolution of 
the ISM represented by the sticky-particle method may 
lead different from that by the SPH method. The SPH 
simulations h ave also still suffere d from the overcooling 
problem (e.g. iKatz fc GunnI Il99il: iNavarro fc Berij | l991 : 



19951 . 



1996 



Navarro fc White! Il994l: iThoul fc Weinberg 
Okamoto et al.ll2005l " ISawala et al.ll201ll ). Degree of dissipa- 
tion, i.e. richness of gas, in merging clumps may also affect 
the properties of the clump-origin bulge. Moreover, struc- 
ture of remnants se ems to depend on the equation of state 
and cooling of gas (|Bournaud et aLlbOllh . Hence, the con- 
tr adiction of our re s ults o f the clump-origin bulge to that 
of lElmegreen et al] (|2008l ) may stem from a difference of 
gaseous physics in the clump merger. The degree of dissipa- 
tion in the clump mergers may depend on initial conditions 
of a simulation model and/or the adopted numerical scheme 
for hydrodynamics (the sticky-particle or SPH code). 

As noted above, the clump-origin bulge in our simu- 
lation has nearly twice larger rotation i n terms of VInax/o"o 



than the re sult oflElmegreen et al] (|200 8'). Rec ently, t h e sim - 
ulations of lAumer et all (120101 " and ICeverino et alj (|201ll 



have also demonstrated that the clumps in clump clusters 
are dynamically rotation supported systems during the mi- 
gration toward the galactic centre althoug h they paid lit- 
tle att ention to resultant bulges. Especially, ICeverino et alj 
(|201ll ) mentioned that clumps which have experienced 
merger with other clumps tend to be highly supported by ro- 
tation and indicated that more massive clumps have larger 
rotation (their §4.2). Although [Ceverino et all (|201ll ) dis- 
cussed kinematics of gas only, their results imply that the 
clump merger can spin up the remnant clumps. Therefore, 
we could expect that the clump-origin bulges in their simu- 
lations would have a significant rotation because the clump- 
origin bulges are the final remnant of the clump mergers. 

Clump-origin bulges form through ^mergers of stellar- 
dumps^ in a disc, whereas classical bulges form through 
^mergers of galaxies'. Formation of the clumps are caused 
by instability of gas in rotating disc, therefore, these clumps 
retain a memory of their discy origin and a s pin in the same 
direction as the disc rotation as shown by ICeverino et alj 
(|201ll ). The bulge formation by the clump mergers takes 
place in the disc plane after orbiting around the galactic 
centre in the same direction, i.e. ^ prograde-prograde merger'. 
As a result, the clump-origin bulges which are remnant of 
the clump merger retain the orbital angular momentum of 
the merger and the spin angular momentum of the progeni- 
tor clumps in their kinematics. Then, the bulges should ac- 
quire the significant rotation of which property is different 
from that of classical bulges although both of them are built 
through merger-origin formation scenarios. 

Our result indicated the nearly exponential surface den- 
sity profile with the Sersic index of n = 1.18. We speculate 
that this exponential profile may stem from the elongated 
bar- like shape of the bulge in our simulation (Fig. |4]). It has 
been observed that barred structu res in spiral galaxies have 
an exponential density profile (e.g. lElmegreen fc ElmegreerJ 
Ohta et alj|l990l ). This was demonstrated also using 



an ex p' 
[l98i;l( 



numerical simulations bv lCombes fc Elmegreeii! (|l993l ) and 
iNoguchil (|l996l fl. 



5 SUMMARY 

We summarize the natures of the clump-origin bulge in our 
simulation: 

• nearly exponential surface density profile, 

• elongated bar-like structure, and boxy shape, 

• significant net rotation, but not cylindrical, 

• old stellar population, 

• high metallicity of stars with vertical gradient and 

• inactive star formation (if there is no late accretion). 

If we compare these natures with t hose observed in typical 
classical bulge and pseudobulges (jKormendv fc KennicuttI 
|2004 ). the former three natures, exponential profile, boxy 
shape and a significant rotation, indicate pseudobulge-like 
signatures. On the other hand, the latter three, old, metal- 
rich stars and weak star formation, represent classical bulge- 
like signatures. Therefore, the clump-origin bulge, which we 
obtained here, can not be simply classified into classical 
bulges nor pseudobulges, but should have unique properties 
as shown in t his paper. 

However, lElmegreen et al] (|2008l ') has suggested with 
numerical simulations that clump-origin bulges should be a 
classical bulge. This inconsistency would be due to adopting 
different initial conditions and/or numerical schemes. As we 
discussed in §3.1.1 and 4, initial setup in simulations would 
control disc diameters, gas density in the disc and the size 
of clump-origin bulges finally, and may affect the natures of 
the clump-origin bulge. 

Hence, our results open up a differ ent possibility of 
clump -origin bulges than the result of lElmegreen et al] 
(|2008l ). In our simulation, the surface density of the 
disc is lower and th e resultant bulge is smaller than 
lElmegreen et al] (|2008l ). Therefore, it is expected that low 
surface brightness clump clusters prefer a clump-origin pseu- 
dobulge and enable the galaxy to avoid the clump-origin 
classical bulge formation. The clump-origin pseudobulge is 
more applicable to the MW bulge which is a pseudobulge- 
like but old metal-rich bulge. The low surface brightness 
clump clusters might had been present but mis sed in the 
observations and lElmegreen fc ElmegreerJ (|2005l ) discussed 
that clump clusters in their observation have a too high sur- 
face density to be the MW progenitor. Thus, to investigate 
the clump-origin bulge formation and its natures may take 
an important role to seek the formation history of the MW. 
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